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The central quantity of density functional theory is the so-called exchange-correlation functional. 
This quantity encompasses all non-trivial many-body effects of the ground-state and has to be ap- 
proximated in any practical application of the theory. For the past 50 years, hundreds of such 
approximations have appeared, with many successfully persisting in the electronic structure com- 
munity and literature. Here, we present a library that contains routines to evaluate many of these 
functionals (around 180) and their derivatives. 
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I. INTRODUCTION 

Density functional theory (DFT) is perhaps one of the most successful theories in Physics and in Chemistry of 
the last half-century [IrS|. It is currently used to predict the structure and the properties of atoms, molecules, and 
solids; it is a key ingredient of the new field of Materials Design, where one tries to create new materials with specific 
properties; it is making its way in Biology as an important tool in the investigation of proteins, DNA, etc. These are 
only a few examples of a discipline that even now, almost 50 years after its birth, is growing at an exponential rate. 

Almost all applications of DFT are performed within the so-called Kohn-Sham scheme 5] , that uses a non-interacting 
electronic system to calculate the density of the interacting system [6[. The Kohn-Sham scheme leeds to the following 
equations. (Hartree atomic units are used throughout the paper, i.e. e 2 = fi = m e = 1.) 



V 2 
— jT" + u cxt(r) + UHartree[n](r) + v xc [n](r) 



ipi(r) =£iipi(r), 



(1) 



where the first term represents the kinetic energy of the electrons; the second is the external potential usually generated 
by a set of Coulombic point charges (sometimes described by pscudopotentials); the third term is the Hartree potential 
that describes the classical electrostatic repulsion between the electrons, 



WHartrccM(r) = Id 
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and the exchange-correlation (xc) potential u xc M is defined by 

SE xc \n] 
v xc [n\(r) 



Sn(r) 



(2) 



(3) 



£>x C [n] is the xc energy functional. Note that by [n] we denote that the quantity is a functional of the electronic 
density, 



t (r)=^|^(r)| 2 , 



(4) 



where the sum runs over the occupied states. 
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FIG. 1: Jacob's ladder of density functional approximations for the xc energy. 



The central quantity of this scheme is the xc energy -ExcM that describes all non-trivial many-body effects. Clearly, 
the exact form of this quantity is unknown and it must be approximated in any practical application of DFT. We 
emphasize that the precision of any DFT calculation depends solely on the form of this quantity, as this is the only 
real approximation in DFT (neglecting numerical approximations that are normally controllable). 

The first approximation to the exchange energy, the local density approximation (LDA), was already proposed by 
Kohn and Sham in the same paper where they described their Kohn-Sham scheme [fj. It states that the value of xc 
energy density at any point in space is simply given by the xc energy density of a homogeneous electron gas (HEG) 
with electronic density n(r). Mathematically, this is written as 



E^P A = /d 3 rn(r)e™ G (n(r)), 



(5) 



where e™ (n) is the xc energy per electron of the HEG. Note that this quantity is a function of n. While the exchange 
contribution e^ EG (n) can be easily calculated analytically, the correlation contribution is usually taken from Quantum 
Monte-Carlo simulations [3,Q|. As defined, the LDA is unique, but, as we will see in the|21 even such precise definition 
can give rise to many different parameterizations. 

During the past 50 years, hundreds of different forms appeared [9j and they are usually arranged in families, which 
have names such as generalized-gradient approximations (GGAs), meta-GGAs, hybrid functionals, etc. In 2001, John 
Perdew came up with a beautiful idea on how to illustrate these families and their relationship [10]. He ordered these 
families as rungs in a ladder that leads to the heaven of "chemical accuracy", and that he christened the Jacob's 
ladder of density functional approximations for the xc energy (see Fig. [I}. Every rung adds a dependency on another 
quantity, thereby increasing the precision of the functional but also increasing the numerical complexity and the 
computational time. 

At the bottom of the ladder we find the LDA, a functional that depends locally on the density only. The second 
rung is occupied by the GGA 



£ - = /^n W e- 



(n(r),Vn(r)). 



(6) 



As one can see, one now adds the gradient of the density, a semi-local quantity that depends on an infinitesimal 
region around r, as a parameter to the energy density. Note that there is a considerable amount of craftsmanship 
and physical/chemical intuition going into the creation of the function e G c (n, Vn), but also a fair quantity of 
arbitrariness. It is therefore not surprising that many different forms were proposed over the years. The same is true 
for the functionals on the next rung, the meta-GGAs 



E K f GA = fd 3 rn(r)e 
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(n(r),Vn(r),V 2 n(r),r(r)). 



(7) 



This time one adds the Laplacian of the density V 2 n(r) and also (twice) the kinetic energy density 

occ 

r(r)=V|V^(r)| 2 . 



(8) 



Note that the meta-GGAs are effectively orbital functionals due to the dependence in r(r). 

The forth rung is occupied by functionals that include the exact-exchange (EXX) contribution to the energy 
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(9) 



These functionals can include the whole £ , ^ xx or only a fraction of it, and E^ xx can be evaluated with the bare 
Coulomb interaction or with a screened version of it. Furthermore, one can add a (semi-)local xc term or one that 
depends on the orbitals. In any case, the fourth rung only includes functionals that depend on the occupied orbitals 
only. Notorious examples of functionals on this rung are the hybrid functionals 

^ mGGA [n(r), Vn(r), V 2 n(r), r(r)] . 
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(10) 



Note that the local part of this functional can be a meta-GGA, a GGA, or even an LDA. 

Finally, on the last rung of Jacob's ladder one finds functionals that depend on the empty (virtual) Kohn-Sham 
orbitals. Perhaps the best known example of these functionals is the random-phase approximation (RPA). 

In a practical DFT calculation one needs typically to evaluate both E xc and v xc . Furthermore, to obtain response 
properties higher derivatives of E xc are required. For example, in first order one can get the electric polarizability, 
the magnetic susceptibility, phonon frequencies, etc., and these usually require the knowledge of the xc kernel 
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Sn(r)5n(r') 
In second order one can obtain, e.g., hyperpolarizabilities, Raman tensors, etc., but these calculations usually require 



fc xc (r,r',r") 
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(12) 



The calculation of these derivatives is fairly straightforward by using basic functional analysis and the chain rule for 
functional derivatives. We give here as an example the case of the xc potential for a GGA 
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Integrating by parts the third term, one finally arrives at 
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Clearly, higher functional derivatives of _E XC involve higher partial derivatives of e xc . 

By now it is clear what a code needs to implement for the functionals of the first three rungs: given n(r), and 
possibly Vn(r), V 2 n(r), and r(r), one needs the xc energy density e xc (r) and all relevant partial derivatives of this 
quantity. This is indeed the information that is provided by Libxc. For hybrid functionals the library also returns, 
besides the semi-local part, also the mixing coefficient a x . Unfortunately, EXX and other functionals of the forth 
and fifth rang are too dependent on the actual numerical representation of the wave-functions and can not be easily 
included in a generic library. 

Finally, we would like to mention that some applications of DFT do not use the Kohn-Sham scheme, but try to 
approximate the kinetic energy functional directly in terms of the density. This approach follows the path that was 
laid down in the 1927 by Thomas and Fermi [lllll2|. and is sometimes referred to as "orbital-free DFT" [13]. In this 
case, we also need a functional form (either an LDA or a GGA) for the kinetic energy, and many of these are present 
in LlBXC . A good review of functionals for the kinetic energy density can be found in Ref. [14J ] . 

Note that even if the discussion above was restricted to spin-compensated systems for the sake of simplicity, all 
functionals of Libxc can also be used with spin-polarization. 



II. SOME HISTORY 

LlBXC started as a spin-off project during the initial development of the (time-dependent) DFT code OCTOPUS fla - 
[l8[ • At that point it became clear that the task of evaluation of the xc functional was completely independent of the 
main structure of OCTOPUS, and could therefore be transformed into a library. The first steps into the development 
of Libxc were taken in September 2006, and the first usable version of the library included a few of the most popular 
LDA and GGA functionals. 

At the same time, the European Theoretical Spectroscopy Facility was carrying a coordinate effort to improve code 
interoperability and re-usability of its software suite, so efforts were almost immediately made to interface those codes 
with Libxc. This catalyzed the development of the library and accelerated its dissemination in the community. 

During the following years, the library was expanded following two main lines: 

(i) To include as many functionals as possible (by now we include around 180 — see IA"|) . This was done for 
several reasons. First, we implemented nearly all of the "old" functionals that played an important role in the 
development of DFT. From this perspective, we can look at Libxc as a "living museum" of the history of this 
important discipline. It also allows users to reproduce old results with little effort. Secondly, many of the functionals 
that are proposed nowadays are implemented in Libxc within weeks from the moment they are published. In this 
way new developments are very rapidly available in several different codes, allowing for these new functionals to be 
quickly tested and benchmarked. Note that we have an agnostic policy, i.e. we try to include the maximum possible 
number of functionals, without making any judgment of value concerning their beauty, elegance, or usefulness. This 
judgment is left for the final user to perform. 

(ii) To include derivatives of the xc energy up to high orders. As we stated before, to perform a standard Kohn- 
Sham calculation one only requires the xc energy functional and its first derivatives. However, higher derivatives are 
essential in order to obtain response properties. Of course, these higher derivatives can be calculated numerically from 
lower-order derivatives, but this procedure tends to introduce unnecessary errors and instabilities in the calculations. 
We therefore implemented derivatives up to the third order for the LDAs and up to second order for the other 
functionals. We note that these derivatives are hand-coded, and not automatically generated from the output of 
symbolic manipulation software |19l |20j . Even if the latter approach is excellent in order to test implementations, 
the automatically generated code is often extremely verbose, inefficient, and unreadable. Note that Libxc includes 
automatic procedures to check the implementation of the analytic derivatives. 

At the beginning, Libxc was used exclusively in the code OCTOPUS. However, since then several other codes from 
both the Solid-State Physics and Quantum Chemistry communities started to use this library. The list of codes that 
use Libxc at the time of writing is as follows (in alphabetical order): 



• Abinit [2lM23j (http://www.abinit.org/) — This is a general purpose plane- wave code. Besides the basic 
functionality, Abinit also includes options to optimize the geometry, to perform molecular dynamics simulations, 
or to calculate dynamical matrices, Born effective charges, dielectric tensors, and many more properties. 



• 



APE [2J] (http://www.tddft.org/programs/APE) — The atomic pseudopotential engine (APE) is a tool for 



generating atomic pseudopotentials within DFT. It is distributed under the GPL and it produces pseudopotential 
files suitable for use with several codes. 

• AtomPAW 25] (http://www.wfu.edu/~natalie/papers/pwpaw/man.html) — The computer program Atom- 
Paw generates projector and basis functions which are needed for performing electronic structure calculations 
based on the projector augmented wave (PAW) method. The program is applicable to materials throughout the 
periodic table. 

• Atomistix ToolKit [H, [27} (http://quantumwise.com) -- This is a software package that uses non- 
equilibrium Green's functions simulations to study transport properties like I-V characteristics of nanoelectronic 
devices. It uses a powerful combination of DFT, semi-empirical tight-binding, and classical potentials. 

• BigDFT [H, HI] ( |http : //inac . cea . f r/L_Sim /BigDFT/ ) — BigDFT is a DFT massively parallel electronic 
structure code using a wavelet basis set. Wavelets form a real space basis set distributed on an adaptive mesh. 
Thanks to its Poisson solver based on a Green's function formalism, periodic systems, surfaces and isolated 
systems can be simulated with the proper boundary conditions. 

• DP [20] (http://www.dp-code.org/) — Dielectric properties (DP) is a Linear Response TDDFT code, in 
frequency-reciprocal and frequency-real space, that uses a plane-wave basis set. 

• Elk (http://elk.sourceforge.net/) — An all-electron full-potential linearized augmented-plane wave (FP- 
LAPW) code with many advanced features. This code is designed to be as simple as possible so that new 



developments in the field of density functional theory (DFT) can be added quickly and reliably. The code is 
freely available under the GNU General Public License. 

ERKALE (http://erkale.googlecode.com) — ERKALE is a quantum chemistry program developed by 
J. Lehtola used to solve the electronic structure of atoms, molecules and molecular clusters. The main use 
of ERKALE is the computation of X-ray properties, such as ground-state electron momentum densities and 
Compton profiles, and core (X-ray absorption and X-ray Raman scattering) and valence electron excitation 
spectra of atoms and molecules. 



• Exciting [HI, HI] ( |http : //exciting- code . org/[ ) — Exciting is a full-potential all-electron DFT package 



• 



based on the linearized augmented plane-wave (LAPW) method. It can be applied to all kinds of materials, 
irrespective of the atomic species involved, and also allows for the investigation of the atomic-core region. 

GPAW [Hill (http sTTTwiki . f ysik . dtu . dk/gpaw) — GPAW is a DFT Python code based on the projector- 
augmented wave (PAW) method and the atomic simulation environment (ASE). It uses real-space uniform grids 
and multigrid methods or atom-centered basis-functions. 



• Hippo [35H37I ] - — This is an electronic structure code, developed by N. Lathiotakis, that implements Reduced 
Density Matrix Functionals for atomic and molecular systems using Gaussian-type orbitals . 

• OCTOPUS |15l - ll8| (http://www.tddft.org/programs/octopus/) — Octopus is a scientific program aimed 
at the ab initio virtual experimentation on a hopefully ever-increasing range of system types. Electrons 
are described quantum-mechanically within density- functional theory (DFT), and in its time-dependent form 
(TDDFT) when doing simulations in time. Nuclei are described classically as point particles. Electron-nucleus 
interaction is described within the pseudopotential approximation. 



• Yambo 38[ ( |http : //www . yambo-code . org/p — Yambo is a Fortran/C code for many-body calculations in 
solid state and molecular physics. The code was originally developed in the Condensed Matter Theoretical 
Group of the Physics Department at the University of Rome "Tor Vergata" by A. Marini. Previous to its 
release under the GPL license, Yambo was known as SELF. 

This diversity of codes is extremely important, because in this way a certain functional can be tested and used 
with a variety of methods in different physical situations. For example, a meta-GGA developed within a certain code 
to get good band-gaps of solids can be immediately used in a different code (often after a simple recompilation) to 
obtain the ionization potential of molecules. 

Libxc is freely available from http://www.tddft.org/programs/Libxc, and it is distributed under the GNU 
Lesser General Public License v3.0. This license not only allows everyone to read, modify, and distribute the code, 
but also allows LlBXC to be linked from close-source codes. The reader is also referred to the web site to obtain more 
information, updated documentation, examples, new versions, etc. 

As most open source projects, we strongly encourage contributions from researchers willing to contribute with the 
implementation of new functionals, higher derivatives, bug corrections, or even bug reports. 

III. AN EXAMPLE 

A. Calling Libxc 

Probably the best way to explain the usage of LlBXC is through an example. The following small program calculates 
the xc energy for a given functional for several values of the density; the available C bindings can be found in header 
file xc.h. More information and examples can be found in the manual and in the header files. 

#include <stdio.h> 
#include <xc.h> 

int mainO 
{ 

xc_func_type func; 

double rho[5] = {0.1, 0.2, 0.3, 0.4, 0.5}; 

double sigma [5] = {0.2, 0.3, 0.4, 0.5, 0.6}; 

double ek[5] ; 



int i, func_id = 1; 

/* initialize the functional */ 

if (xc_func_init(&func, func_id, XC_UNPOLARIZED) != 0) { 

fprintf (stderr , "Functional '7,d' not f ound\n" , func_id) ; 

return 1 ; 
} 

/* evaluate the functional */ 
switch(f unc . inf o->f amily) 
{ 
case XC_FAMILY_LDA: 

xc_lda_exc(&f unc, 5, rho, ek) ; 

break; 
case XC_FAMILY_GGA: 
case XC_FAMILY_HYB_GGA: 

xc_gga_exc(&f unc, 5, rho, sigma, ek) ; 

break; 
} 

for(i=0; i<5; i++) { 

printf ("°/.lf 7.1f\n", rho [i] , ek[i]); 
} 

/* free the functional */ 
xc_func_end(&func) ; 
} 

The functionals are divided in families (LDA, GGA, etc.). Given a functional identifier, f unc_id, the functional is 
initialized by xc_func_init, and evaluated by xc_XXX_exc, which returns the energy per unit volume (ek). Finally, 
the function xc_f unc_end cleans up. We note that we follow the convention used in Quantum Chemistry and, instead 
of passing the full gradient of the density, we use the variable 

cr(r) = Vn(r) • Vn(r) . (15) 

Converting between partial derivatives with respect to a and Vn is trivially done with the use of the chain rule. All 
the quantities passed to and returned by the library are in atomic units. 

Fortran 90 bindings are also included in Libxc. These can be found in the file libxejnaster .F90. In general, 
calling Libxc from Fortran is as simple as from C. Here is the previous example in Fortran: 

program lxctest 

use xc_f 90_types_m 
use xc_f 90_lib_m 

implicit none 

TYPE(xc_f90_pointer_t) : : xc_func 

TYPE(xc_f 90_pointer_t) : : xc_info 

real (8) 

real (8) 

real (8) : : ek(5) 

integer :: i, func_id 



rho(5) = (/0.1, 0.2, 0.3, 0.4, 0.5/) 
sigma(5) = (/0.2, 0.3, 0.4, 0.5, 0.6/) 



func_id = 1 

! initialize the functional 

call xc_f90_func_init(xc_func, xc_info, func_id, XC_UNP0LARIZED) 

! evaluate the functional 



select case (xc_f 90_inf o_f amily(xc_inf o)) 
case (XC_FAMILY_LDA) 

call xc_f 90_lda_exc(xc_f unc, 5, rho(l), ek(l)) 
case (XC_FAMILY_GGA , XC_FAMILY_HYB_GGA) 

call xc_f 90_gga_exc(xc_f unc, 5, rho(l), sigma(l) , ek(l)) 
end select 

do i = 1, 5 

write(*,"(F8.6,lX,F9.6)") rho(i), ek(i) 
end do 

! free the functional 

call xc_f 90_func_end(xc_func) 

end program lxctest 



B. The info structure 

Besides the mathematical formulas necessary to evaluate the functional and its derivatives, Libxc includes a 
considerable amount of metadata that is quite useful for both the calling program and for the end user. This 
information is contained for each functional in the structure xc_f unc_inf o_type. The relevant part of this structure 
for the end user is defined as 



/* flags that can be used in info. 
#define XC_FLAGS_HAVE_EXC 
#define XC_FLAGS_HAVE_VXC 
#define XC_FLAGS_HAVE_FXC 
#define XC_FLAGS_HAVE_KXC 
#define XC_FLAGS_HAVE_LXC 
#define XC_FLAGS_1D 
#define XC_FLAGS_2D 
#define XC_FLAGS_3D 
#define XC_FLAGS_STABLE 
#define XC FLAGS DEVELOPMENT 



flags */ 



(1 


« 


0) 


/* 


1 


*/ 


(1 


« 


1) 


/* 


2 


*/ 


(1 


« 


2) 


/* 


4 


*/ 


(1 


« 


3) 


/* 


8 


*/ 


(1 


« 


4) 


/* 


16 


*/ 


(1 


« 


5) 


/* 


32 


*/ 


(1 


« 


6) 


/* 


64 


*/ 


(1 


« 


7) 


/* 


128 


*/ 


(1 


« 


9) 


/* 


512 


*/ 


(1 


« 


10) 


/* 


1024 


*/ 



typedef struct{ 
int number ; 
int kind; 



/* indentifier number */ 

/* XC_EXCHANGE and/or XC_C0RRELATI0N */ 



char *name ; 
int family; 
char *refs; 



/* name of the functional, e.g. "PBE" */ 

/* type of the functional, e.g. XC_FAMILY_GGA */ 

/* references */ 



int flags; /* see above for a list of possible flags */ 



} xc_f unc_inf o_type; 

For example, for the Slater exchange functional, this structure is defined as 

const XC(func_inf o_type) XC(func_inf o_lda_x) = { 
XC_LDA_X, 
XC_EXCHANGE, 
"Slater exchange", 
XC_FAMILY_LDA, 
"PAM Dirac, Proceedings of the 

Cambridge Philosophical Society 26, 376 (1930) \n" 
"F Bloch, Zeitschrift fuer Physik 57, 545 (1929)", 



XC_FLAGS_3D I XC_FLAGS_HAVE_EXC I XC_FLAGS_HAVE_VXC I XC_FLAGS_HAVE_FXC 
I XC_FLAGS_HAVE_KXC, 

>; 

Note that the references are separated by a newline. The user of the library can access the information in the following 
way: 

#include <stdio.h> 
#include <xc.h> 

int mainO 
{ 

xc_func_type func; 

xc_func_init(&func, XC_GGA_X_B88, XC_UNPDLARIZED) ; 

printfC'The functional '°/ s' is defined in the ref erence(s) : \n°/s\n" , 
func . inf o->name , func . inf o->ref s) ; 

xc_f unc_end(&f unc) ; 



IV. CONCLUSIONS AND THE FUTURE 

LlBXC is by now seven years old, and the code is quite stable and in use by hundreds of scientists around the 
world. We are committed to continue the development of the library in the future, mainly by following the two 
lines mentioned before: include all functionals, and their derivatives of the highest possible order. By now essentially 
all LDA and GGA functionals ever proposed in the literature are already included in the library. Unfortunately, 
important gaps still remain especially in the meta-GGAs and hybrid functionals. We will fill in those gaps in the near 
future. It is not clear if Libxc will ever include all functionals ever developed (especially, as several new functionals 
come out every year) , but we will of course try. . . 
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Appendix A: Available functionals 

For convenience, functionals in Libxc are divided in exchange, correlation, and exchange-correlation functionals. 
This division is sometimes arbitrary, but it is often useful for the end user (who may fancy strange mixtures of 
functionals) and for the code developers. 

Below we give a complete list of the functionals the library currently knows about. We note that many of these 
were tested against reference implementations, when available, or against published results. The functionals have a 
label that is used to identify the functional inside Libxc. Furthermore, the third column indicates the time in seconds 
required for 50,000,000 evaluations of the (spin-unpolarized) xc potential in a single core of an Intel Core 2 processor 
running at 2.83 GHz. As one can see, the typical execution time within Libxc is much smaller than the total time 
required in the whole calculation. 

TABLE I: Functionals available in Libxc. 
LDA Functionals 



LDA Exchange 

XCJLDAJt 

XC_LDA_K_2D 

XC_LDA_X_1D 

LDA Correlation 

XC_LDA_C.WIGNER 

XC_LDA_C_RPA 

XCXDA_C_HL 

XCXDA.C.GL 

XCJLDA_C_XALPHA 

XCJLDA.C.VWN 

XCXDA_C.VWN_RPA 

XC_LDA.C_PZ 

XC_LDA_C.OB_PZ 

XCJLDA_C_PW 

XCJLDA_C_PW_RPA 

XCXDA_C.OB_PW 

XCXDA_C_2D_AMGB 

XCJLDA_C_2D_PRM 

XC_LDA_C_vBH 
XCXDA_C_1D_CSC 

XCJLDA.CJVILI 

XCXDA_C.ML2 

XC_LDA_C_GOMBAS 



LDA exchange 
Slater exchange in 2D 
Slater exchange in ID 



Wigner parametrization 
Random Phase Approximation 
Hedin & Lundqvist 
Gunnarsson & Lundqvist 
Slater's Xa (X-alpha) 
Vosko, Wilk, & Nussair 
Vosko, Wilk, & Nussair (RPA) 
Perdew & Zunger 

Ortiz & Ballone (PZ parametrization) 
Perdew & Wang 

Perdew & Wang fit to the RPA energy 
Ortiz & Ballone (PW parametrization) 
Attaccalite, Moroni, Gori-Giorgi, and 
Bachelet (LDA for 2D systems) 
Pittalis, Rasanen, and Marques (LDA for 
2D systems) 
von Barth & Hedin 

Casula, Sorella & Senatore (LDA correla- 
tion for ID systems) 

Modified LDA (version 1) of Proynov and 
Salahub 

Modified LDA (version 2) of Proynov and 
Salahub 
Gombas 



6.92 


[40, 41] 


9.36 




3382 


f42] 


6.55 


[43] 


11.96 


[44] 


9.86 


[45] 


12.19 


[46] 


9.33 




24.14 


[47] 


18.21 


[47] 


8.90 


[48] 


12.08 


[8,48,49] 


17.63 


[50] 


28.16 


[50] 


12.43 


[8,49,50] 


10.44 


[51] 


11.35 


[52] 


10.40 


[53] 


28.85 


[54] 


37.02 


[55] 


51.99 


[55] 


10.81 


[56] 



LDA Exchange-Correlation 

XC_LDA_XC_TETER93 Teter 1993 



9.22 [57] 



LDA Kinetic Energy 

XCJLDAJK.TF 
XCXDA_K_LP 



Thomas- Fermi kinetic energy 6.74 [11, 12] 

Lee and Parr Gaussian ansatz for the ki- 7.99 [58] 
netic energy 



GGA Functionals 



GGA Exchange 

XC.GGA_X_PBE 
XC.GGA_X_PBE_R 

XC.GGA_X_MPBE 
XC.GGAJCJCPBE 
XC_GGA_X_B86 
XC.GGA_X_B86.MGC 

XC.GGA_X_B88 

XC.GGA_X_G96 

XC_GGA_X_PW86 

XC_GGA_X_PW91 

XC.GGA_X_OPTX 

XC.GGA_X_DK87.R1 

XC.GGA_X_DK87.R2 

XC_GGA_X_LG93 

XC.GGA_X_FT97.A 

XC.GGA_X_FT97.B 

XC.GGA_X_PBE_SOL 

XC.GGA_X_RPBE 
XC.GGA_X_WC 



Perdew, Burke & Ernzerhof exchange 


13.61 


[59,60] 


Perdew, Burke & Ernzerhof exchange 


19.27 


w 


(revised) 






Adamo & Barone modification to PBE 


19.39 


[62] 


Extended PBE by Xu & Goddard 111 


18.90 


[63] 


Becke 86 Xalfa,beta,gamma 


18.64 


[64,65] 


Becke 86 Xalfa,beta,gamma (with mod. 


24.86 


[64,66] 


grad. correction) 






Becke 88 


20.97 


[67] 


Gill 96 


15.57 


[68] 


Perdew & Wang 86 


32.31 


[69] 


Perdew & Wang 91 


48.02 


[70] 


Handy & Cohen OPTX 01 


18.21 


[71] 


dePristo & Kress 87 (version Rl) 


19.10 


[72] 


dePristo & Kress 87 (version R2) 


24.11 


[72] 


Lacks & Gordon 93 


32.74 


[73] 


Filatov & Thiel 97 (version A) 


22.05 


[74] 


Filatov & Thiel 97 (version B) 


26.86 


[74] 


Perdew, Burke & Ernzerhof exchange (for 


19.33 


[75] 


solids) 






Hammer, Hansen & Norskov (PBE-like) 


17.14 


[76] 


Wu & Cohen 


17.91 


[77] 
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XC_GGA_X_AM05 
XC.GGA_X_PBEA 
XC.GGA_X_mPW91 
XC.GGA_X_2D_B86.MGC 

XC.GGA_X_BAYESIAN 

XC_GGA_X_PBE_JSJR 

XC.GGA_X_OPTB88.VDW 

XC.GGA_X_PBEK1_VDW 

XC.GGA_X_OPTPBE_VDW 

XC.GGA_X_RGE2 

XC_GGA_X_RPW86 

XC.GGA_X_KT1 

XC.GGA_X_HERMAN 

XC.GGA_X_LBM 

XC.GGA_X_OL2 

XC.GGA_X_MB88 
XC.GGA_X_APBE 

XC.GGA_X_HTBS 
XC.GGA_X_AIRY 
XC.GGA_X_LAG 
XC.GGA_X_SOGGAll 

XC.GGA_X_C09X 



Armiento & Mattsson 05 exchange 
Madsen 07 

mPW91 of Adamo & Barone 
Becke 86 with modified gradient correction 
for 2D 

Bayesian best fit for the enhancement 
factor 

Reparametrized PBE by Pedroza, Silva & 
Capelle 

opt-Becke 88 for vdW 
Reparametrized PBE for vdW 
Reparametrized PBE for vdW 
Regularized PBE 
Refitted Perdew & Wang 86 
Keal and Tozer, version 1 
Herman Xalphabeta GGA 
van Leeuwen & Baerends modified 
Exchange form based on Ou-Yang and 
Levy v.2 

Modified Becke 88 for proton transfer 
mu fixed from the semiclassical neutral 
atom 

Haas, Tran, Blaha, and Schwarz 
Constantin et al based on the Airy gas 
Local Airy Gas 

Second-order generalized gradient approx- 
imation 2011 

C09x to be used with the VdW of Rutgers- 
Chalmers 



61.38 


[78,79] 


30.98 


f80] 


47.19 


[81] 


13.49 


[82] 


18.91 


[83] 


18.71 


[84] 


25.24 


[85] 


18.72 


[85] 


40.62 


[85] 


14.65 


[86] 


30.61 


[87] 


23.96 


[88] 


17.44 


[89,90] 


30.98 


[211 


15.69 


[92,93] 


24.86 


[94] 


18.71 


[95] 


26.17 


[96] 


47.48 


[97] 


49.75 


[98] 


20.40 


[991 



19.; 



GGA Correlation 
XC.GGA_C.PBE 
XC_GGA_C_XPBE 
XC.GGA.CJLYP 
XC_GGA_C_P86 
XC_GGA_C_PBE_SOL 

XC.GGA_C.PW91 
XC.GGA_C_AM05 
XC_GGA_C_LM 
XC_GGA_C_PBE_JRGX 

XC.GGA_C_RGE2 

XC.GGA_C.WL 

XC.GGA_C.WI 

XC.GGA_C.WI0 

XC_GGA_C_APBE 

XC_GGA_C_SOGGAll 

XC_GGA_C_SOGGAll_X 



Perdew, Burke & Ernzerhof correlation 39.69 

Extended PBE by Xu & Goddard III 43.85 

Lee, Yang & Parr 25.88 

Perdew 86 51.36 

Perdew, Burke & Ernzerhof correlation 55.32 

SOL 

Perdew & Wang 91 

Armiento & Mattsson 05 correlation 

Langreth & Mehl 

Reparametrized PBE by Pedroza, Silva . 

Capelle 

Regularized PBE 40.55 

Wilson & Levy 20.95 

Wilson & Ivanov 16.48 

Wilson & Ivanov initial version 19.42 

mu fixed from the semiclassical neutral 40.86 

atom 

Second-order generalized gradient approx- 40.67 [991 

imation 2011 

to be used with 37.05 [1081 

XC_HYB_GGA_X_SOGGAll_X 



[59,60] 
[63] 

[101, .102] 
[103J 

[751 



51.17 [70. 1041 

25.88 [78] ~ 

42.87 [105] 

55.74 [841 



[86] 

[106] 

[107] 

[107] 
[95] 



GGA Exchange-Correlat 

XC_GGA_XC_LB 

XC_GGA_XC_HCTH_93 

XC_GGA_XC_HCTH_120 

XC_GGA_XC_HCTH_147 

XC_GGA_XC_HCTH_407 

XC_GGA_XC_EDF1 

XC_GGA_XC_XLYP 
XC_GGA_XC_PBE1W 

XC.GG A_XC_MP WLYP 1 W 



ion 

van Leeuwen & Baerends 

HCTH functional fitted to 93 molecules 

HCTH functional fitted to 120 molecules 

HCTH functional fitted to 147 molecules 

HCTH functional fitted to 147 molecules 

Empirical functional from Adamson, Gill, 

and Pople 

XLYP functional 

PBE1W (functional fitted for water) 

mPWLYPlw (functional fitted for water) 



29.27 
155.10 
151.57 
143.57 
155.60 

69.03 



[109] 
[110] 

[111] 
[112] 
[113] 



106.30 [114] 

94.72 [115J 

106.74 [1151 



11 



XC_GGA_XC_PBELYP1W 

XC.GGA_XCJKT2 

XC.GGA_XC.TH_FL 

XC.GGA_XC.TH_FC 

XC.GGA_XC.TH_FCFO 

XC.GGA_XC.TH_FCO 

XC.GGA_XC.TH1 

XC.GGA_XC.TH2 

XC_GGA_XC_TH3 

XC.GGA_XC.TH3 

GGA Kinetic Energy 

XC.GGA.K_VW 

XC.GGA_K.GE2 

XC_GGA_K_GOLDEN 

XC_GGA_K_YT65 

XC_GGA_K_BALTIN 
XC_GGA_K_LIEB 

XC_GGA_K_ABSR1 

XC_GGA_K_ABSR2 

XC_GGA_K_GR 



XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 
XC.GGA. 



K.LUDENA 

K.GP85 

K.LLP 

K_FR_B88 

K_FR_PW86 

K.PEARSON 

K.0L1 

K.OL2 

K.DK 

K.PERDEW 

K.VSK 

K.VJKS 

K.ERNZERHOF 

K.LC94 

K.APBE 



XC.GGA.K.THAKKAR 
XC.GGA_K.TW1 
XC.GGA_K_TW2 
XC_GGA_K_TW3 

XC_GGA_K_TW4 

MetaGGA Functionate 



PBELYP1W (functional fitted for water) 79.12 

Keal and Tozer, version 2 31.73 

Tozer and Handy v. FL 94.64 

Tozer and Handy v. FC 247.35 

Tozer and Handy v. FCFO 372.59 

Tozer and Handy v. FCO 355.16 

Tozer and Handy v. 1 383.44 

Tozer and Handy v. 2 352.52 

Tozer and Handy v. 3 332.24 

Tozer and Handy v. 4 337.86 



[115] 

f88] 

[116] 

[116] 

[116] 

[116] 

[117] 

[118] 

[119] 

[119] 



von Weiszaecker correction to Thomas- 
Fermi 

Second-order gradient expansion of the ki- 
netic energy density 

TF-lambda-vW form by Golden (I = 
13/45) 

TF-lambda-vW form by Yonei and 
Tomishima (/ = 1/5) 
TF-lambda-vW form by Baltin (/ = 5/9) 
TF-lambda-vW form by Lieb (I = 
0.185909191) 

gamma-TFvW form by Acharya et al [g = 
1 - 1A12/N 1/3 ] 

gamma-TFvW form by Acharya et al [g = 
1 - 1.332/_V 1/3 ] 

gamma-TFvW form by Gazquez and 
Robles 

gamma-TFvW form by Ludefia 
gamma-TFvW form by Ghosh and Parr 
Becke 88 

Fuentealba & Reyes (B88 version) 
Fuentealba & Reyes (PW86 version) 
Pearson 1992 
Ou-Yang and Levy v.l 
Ou-Yang and Levy v. 2 
DePristo and Kress 
Perdew 

Vitos, Skriver, and Kollar 
Vitos, Johansson, Kollar, and Skriver 
Ernzerhof 

Lembarki & Chermette 
mu fixed from the semiclassical neutral 
atom 

Thakkar 1992 

Tran and Wesolowski set 1 (Table II) 
Tran and Wesolowski set 2 (Table II) 
Tran and Wesolowski set 3 (Table II) 
Tran and Wesolowski set 4 (Table II) 



18.49 [120] 
16.08 [121 . 1221 
16.39 [123] 
13.29 [1241 



13.24 
14.68 



125] 
"1261 



17.54 [12JJ 
17.01 [127] 
13.89 [1281 



18.72 
18.19 
22.69 
19.97 
32.24 
13.30 
20.43 
20.26 
18.64 
16.21 
17.05 
18.61 
18.72 
46.63 
19.48 



[129] 
[130] 

[121] 

[92] 

[92] 

[132, 

[93] 

[93] 

[134] 

[135] 

[136] 

[137] 

[138] 

[139] 

[95] 



24.92 [140] 

18.70 [141] 

18.30 [141] 

19.01 [141] 

19.97 [1411 



MetaGGA Exchange 

XC_MGGA_X_LTA 

XC_MGGA_X_TPSS 

XC_MGGA_X_TAU_HCTH, 

XC_MGGA_X_GVT4 

XC_MGGA_X_M06L 

XC_MGGA_X_BR89 

XC_MGGA_X_BJ06 

XC_MGGA_X_TB09 



Local tau approximation 

Perdew, Tao, Staroverov & 

exchange 

Tau HCTH 

GVT4 (exchange part of VSXC) 

M06-L 

Becke-Roussel 89 

Becke & Johnson 06 

Tran & Blaha 09 



Scuseria 



24.07 
28.88 

23.45 

26.85 
29.74 
49.39 
35.28 
31.96 



142] 
"143, 



145] 
146] 
"147. 
149] 
150] 
"151] 
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XCJdGGA_X_RPP09 Rasanen, Pittalis & Proetto 09 44.16 [152] 

XCJVIGGA_X^D_PRHG07 Pittalis-Rasanen-Helbig-Gross 2007 54.18 [153] 

XC_MGGA_X^D_PRHG07_PRP10 Corrected Pittalis-Rasanen-Helbig-Gross 53.94 [153] 

2010 



MetaGGA Correlation 

XCJVIGGA_C_TPSS Perdew, Tao, Staroverov & Scuseria 184.60 [143, 144] 

correlation 
XCJVIGGA_C_VSXC VSXC (correlation part) 148.84 [1461 

Hybrid Functionals 



Hybrid Exchange 

XCJHYB_GGA_X_SOGGAll_X 



Hybrid based on SOGGA11 form 
XC_GGA_C_SOGGAll_X) 



(see also 21.13 [1081 



Hybrid Exchange-Correlation 
XCJHYB_GGA_XC_B3PW91 
XC_HYB_GGA_XC.B3LYP 
XC_HYB_GGA_XC_B3P86 
XC_HYB_GGA_XC_03LYP 
XC_HYB_GGA_KC.PBEH 
XC_HYB_GGA_XC_X3LYP 
XC_HYB_GGA_XC_B1WC 
XC_HYB_GGA_XC.B97 
XC_HYB_GGA_XC.B97.1 
XCJHYB_GGA_XC.B97.2 
XCJHYB_GGA_XC.B97.K 
XC_HYB_GGAJCC.B97.3 
XC_HYB_GGA_XC.B1LYP 
XC_HYB_GGA_XC.B1PW91 
XCJHYB_GGA_XC.mPWlPW 
XCJlYB_GGA_XC_mPW3PW 
XC_HYB_GGA_XC.mPW3LYP 
XC_HYB_GGA_XC.SB98.la 
XC_HYB_GGA_XC.SB98.lb 
XC_HYB_GGAJCC.SB98.lc 
XC_HYB_GGA_XC.SB98.2a 
XC_HYB_GGA_XC_SB98_2b 
XC_HYB_GGA_XC.SB98.2c 
XCJHYB_GGAJtC_mPWlK 



The original hybrid proposed by Becke 


99.90 


154] 


The (in)famous B3LYP 


93.55 


155] 


Perdew 86 hybrid similar to B3PW91 


105.17 




hybrid using the optx functional 


78.60 


1561 


PBEH, also known as PBE0 


76.12 


157] 


X3LYP 


120.61 


114] 


B1WC 


78.94 


158] 


Becke 97 


158.49 


65] 


Becke 97-1 


154.97 


110] 


Becke 97-2 


151.12 


159] 


Becke 97-K, Boese-Martin for Kinetics 


158.38 


160] 


Becke 97-3 


155.92 


161] 


B1LYP 


52.18 


162] 


B1PW91 


69.92 


162] 


mPWIPW 


99.96 


81] 


mPW3PW of Adamo & Barone 


110.82 


8JJ 


mPW3LYP 


112.92 


163] 


SB98 (la) 


149.81 


164] 


SB98 (lb) 


156.87 


164] 


SB98 (lc) 


158.72 


164] 


SB98 (2a) 


139.63 


164] 


SB98 (2b) 


148.48 


164] 


SB98 (2c) 


158.78 


164] 


mPWIK 


100.53 


165 



[11 
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